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Pairs of trapped atoms can be associated to make a diatomic molecule using a time dependent 
magnetic field to ramp the energy of a scattering resonance state from above to below the scattering 
threshold. A relatively simple model, parameterized in terms of the background scattering length 
and resonance width and magnetic moment, can be used to predict conversion probabilities from 
atoms to molecules. The model and its Landau-Zener interpretation are described and illustrated 
by specific calculations for 23 Na, 87 Rb, and 133 Cs resonances. The model can be readily adapted to 
Bose-Einstein condensates. Comparison with full many-body calculations for the condensate case 
show that the model is very useful for making simple estimates of molecule conversion efficiencies. 

PACS numbers: 



I. INTRODUCTION 

The production and characterization of cold, trapped 
molecules and quantum degenerate molecular gases is 
an important contemporary challenge in physics. Con- 
siderable progress towards this goal has occurred re- 
cently, using the properties of magnetically tunable Fes- 
hbach resonance states to convert atoms in an ultracold 
Bose 0, H H H or Fermi H 0, S H E3 gas to bound 
state molecules comprised of a pair of atoms. This pa- 
per describes a simple model for understanding the ba- 
sic molecular physics underlying the conversion of atom 
pairs to diatomic molecules by a time-dependent ramp 
of a scattering resonance across the collision threshold. 
We have in mind the case of two atoms confined in an 
individual cell of an optical lattice, as well as a trapped 
quantum degenerate gas of bosons. Although we do not 
treat quantum degenerate Fermi gases here, the 2-body 
physics we discuss is applicable to cither fermionic or 
bosonic atoms. The molecular physics of threshold scat- 
tering of alkali species is now well-characterized fill . We 
will show how the model proposed by Mies et al. |l2l | gives 
a very simple way to estimate molecule production rates 
in terms of only three parameters that characterize the 
molecular physics and scattering. We also show that the 
model makes predictions that are quantitatively similar 
to the results of full many-body calculations for a trapped 
Bose-Einstein condensate(BEC) using the methods of 
Kohler et al. [H Q E Eg . 

Figure 1 illustrates a Feshbach ramp by showing the 
time-dependent energy levels for a case to be discussed 
below using the model of Ref. ^3 ■ The idea is to take 
a quasibound molecular resonance state whose energy e„ 
is well above the threshold energy e = for the two- 
body collision, and sweep the resonance state energy be- 
low threshold by a time-dependent change in its position. 
Here we consider the case where this is done by tuning 
a magnetic field, although it is also possible to use an 



optical frequency electromagnetic field. If the sweep rate 
is sufficiently slow, then the initially populated ground 
state \v = 0} of a trapped atom pair adiabatically evolves 
to an atom pair in the molecular resonance state \n), 
now located below threshold and stable with respect to 
decay into two atoms. For two atoms in a trap, the 
transfer of population to the final states is accurately de- 
scribed by a Landau-Zener model. In general, the molec- 
ular physics of scattering resonances for alkali species 
is now well-understood, with remarkable agreement be- 
tween coupled channels theoretical models and measure- 
ment 0, 0, 0| . The molecular physics associated with 
an isolated resonance is characterized by three key pa- 
rameters: the background scattering length A^g away 
from resonance, the resonance width A„, and the mag- 
netic moment difference between the resonant state \n) 
and the separated atoms proportional to s n = de n /dB. 
These parameters can, in principle, be determined from 
cither calculation or measurement. In the Landau-Zener 
model the molecule formation probability depends only 
on the product At g A n , and not on s n . 

Ref. 01 showed the connection between the conven- 
tional two-body scattering picture of the resonance in- 
teracting with a continuum and the corresponding prop- 
erties involving discrete quantized levels in a trap. The 
paper also showed how to modify the model when the 
atoms are in a Bose-Einstein condensate (BEC). Here we 
will summarize the key results of Ref. [Tj] , following its 
notation, and extend and apply the method to several 
specific cases. Section [H] considers the collision of two 
atoms in free space, influenced by a single tunable reso- 
nance level. Section IIIII shows how to adapt this picture 
to time-dependent dynamics in a trap when the mag- 
netic field ramps the resonance energy across the collision 
threshold. Section IIVI describes the many-body adapta- 
tion needed for resonance ramps in a Bose-Einstein con- 
densate. Section [V] gives results for specific cases, and 
compares the simple model predictions with those of full 
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FIG. 1: Energy levels for a 6-level model (one tunable res- 
onant state and five trap levels, omitting higher ones) of a 
time-dependent Feshbach resonance ramp for two 87 Rb F= 1, 
M= +1 atoms in a wj, = 27r(100 kHz) isotropic harmonic 
trap. Energies are relative to the energy of two separated 
atoms, which is taken to define e = 0. The dashed lines show 
the diabatic (or "bare") energy levels before the interactions 
between the trap levels \v) and bound state \ri) are turned on. 
The solid lines show the adiabatic (or "dressed") eigenener- 
gies of the full Hamiltonian. The diabatic resonance state 
energy e n is ramped linearly in time by changing the mag- 
netic field B at a speed of 0.5 T/s, such that e n = at time 
t = 0. The time scale is in units of u}^ 1 = 1.16/xs. A pair 
of atoms initially in the ground trap state \v = 0) for time 
t <C — tLifr 1 can evolve adiabatically to the bound molecular 
state |n) for t S> uj^ 1 . The model resonance state \n) has 
a width A n = 0.017 mT and a slope (see text) s n /h — 38 
MHz/mT. This case corresponds to an actual 87 Rb resonance 
near 100.76 mT El. 



many-body calculations. 



II. THRESHOLD RESONANT COLLISIONS IN 
FREE SPACE 

Mies et al. E3 showed how the full coupled channels 
quantum scattering problem in the case of a set of closed 
channels with isolated resonance levels interacting with a 
single scattering continuum can be reduced to an effective 
two-channel configuration interaction (CI) problem. The 
wavefunction at collision kinetic energy e in a magnetic 
field of strength B is represented by the following CI 
wavefunction: 

^ n (e,B,R) = \l)F 1 (e,R,B) + \n)<j )n (R)A n (e 1 B). (1) 

Here R is interatomic separation, |1) represents the 
internal spin variables of the open collision chan- 
nel, F\(e, R, B) is the open-channel scattering function, 
\n)4> n (R) represents the resonant closed-channel bound 
state with internal spin state \n) and motional state 
4> n (R), and A n (e,B) is the amplitude of the resonant 



state. The wavefunction is the solution to the Hamil- 
tonian 



H n (B) 



|1)<1|{T- 
\n)(n\{T 

{|i)H + 



-U bg (R)} + 
\-U n (B,R)} + 
n)(l\}W nA (R) 



(2) 



where T is the radial kinetic energy operator. The |1) 
and \n) basis states are "bare" or diabatic states before 
the interchannel coupling W n ,i(R) is taken into account. 
Ref. 01 demonstrated how the potentials Ub g (R) and 
U n (R) and the coupling operator W n ,i(R) can be re- 
placed by convenient effective potentials, or pseudopo- 
tcntials, and an effective coupling operator. These effec- 
tive potentials and operator reproduce all near-threshold 
scattering properties. 

The amplitude A n (e, B) of the resonant state is in- 
dependent of R, and the bound state radial vibrational 
function <fi n (R) is a solution to the single channel diabatic 
bound state eigenvalue problem: 



2n dR 2 



U n {B, R)\<t> n {R) = e n (B)MR) . (3) 



where U n {B,R) is the closed channel effective potential 
and is the reduced mass of the pair of atoms. The func- 
tion 4> n (R) is just an ordinary vibrational wavefunction, 
the shape of which is independent of B. 

The effect of the magnetic field is merely to shift the 
position of the bound state energy e„ (B) . Let B n define 
the resonance field where the diabatic bound state crosses 
the open channel threshold at e = 0, i. e., e n (B n ) = 0. 
Note that we set the open channel threshold energy, that 
is, the energy of the two separated atoms, to be inde- 
pendent of B. This emphasizes that what is important is 
the energy of the resonance state relative to the energy 
of the two separated atoms. Then for B close to B n , 



e n (B) = s n (B-B n ) 



(4) 



where the slope s n = de n (B) / dB is a constant equal 
to the difference in magnetic moments of the resonance 
state and two separated atoms. Assuming that the mag- 
netic moments do not vary over the small range of B 
considered, varying the magnetic field simply causes the 
resonance energy to ramp linearly across threshold while 
\n)<j) n (R) remains unchanged. 

The open channel scattering solution has the asymp- 
totic form at large R: 



Fx{t,R,B) 



2fi 
nTi 2 k 



sm(kR + £(B)) 



where k = \/2[ie/h and 



(5) 



(6) 



where £t> 9 and £ n (-B) are the respective background and 
resonance contributions to the phase shift £. When B is 
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such that the resonance energy is tuned far from thresh- 
old, the effect of the resonance is negligible, and the scat- 
tering reduces to the background scattering by the open 
channel effective potential U bg (R): 



h z d 2 

2n dR 2 



+ U h& {R)](j} e {R) = e<j> e (R) 



(7) 



where e > is the collision kinetic energy. The asymp- 
totic solution of Eq. Q has the form of Eq. © with 
£ = S,bg- A key collision parameter is the background 
scattering length A bg : 



A bg = lim 



fc^o 2k 



(8) 



There is also a resonance contribution to the phase 
shift, which is important when B is tuned so that e n (B) 
is near threshold: 



tan£ n (e, B) 



r„(e) 



2[e ~ e n (B) - 5e n (B)] ' 



(9) 



where T n (e) is the decay width and Se n (B) a shift from 
the "bare" resonance crossing point (see below). The 
decay width is given by the Fermi golden rule: 



r„(e) ee 27r\V n (e)\ 2 = 2^ n |W n , 1 |0 e )| 2 



(10) 



The Wigner threshold law for s-waves requires that 
r„(e) oc k as k — > 0. 

When k is very small, such that kA <C 1, the definition 
of the scattering length A, 

A(B) = lim = A bg + A n (B) . (11) 

using the k — > limit of Eq. @ , gives the common for- 
mula 



where 



A(B) = A bg 1 - 



B — B„ 



r«(e) 



SB, 



2kA 



bg 



(12) 



(13) 



The width A„ is a constant independent of k close to 
threshold, since the ratio T n (e)/k is independent of k 
due to the threshold law. Note that A bg and s n can be 
positive or negative, so that A„ has the same sign as the 
product A bg s n . 

Using the methods of Ref . [2(j , the shift can be shown 
to be 



SB n = 



r(l — r) 
l + (l-r) ; 



■A„, 



(14) 



where r = A bg /A, and A = 0.477989(2^C 6 /^ 2 ) 1 / 4 is the 
mean scattering length defined by Gribakin and Flam- 
baum 21.1. for a collision of two atoms with reduced mass 



/j, in a van der Waals potential with coefficient C§. The 
trap CI model discussed in the next Section implicitly 
includes the shift, if a complete set of states for the back- 
ground potential U bg is included in the basis. It is es- 
pecially important if A bg is large and positive to include 
the last bound state of U bg in the basis, since it can make 
a large contribution to the shift (see Section IV CI below') . 
However, the final populations in the CI model after the 
ramp arc completely independent of the shift, although 
the detailed dynamics during the ramp can be strongly 
modified when the shift is taken into account. For the 
purpose of this paper, which is primarily concerned with 
the final populations, we can safely ignore the shift. 



III. THRESHOLD RESONANT COLLISIONS IN 
A TRAP 



The free-space picture described in Section [H] needs 
to be modified when a collision occurs in a trap, where 
the atoms are no longer free to separate to arbitrarily 
large R. Their motion is confined to the trapping vol- 
ume, and the free scattering states must be replaced by 
the qua ntized eigenstates of the interacting atoms in the 
trap[22, |U |24|. The center of mass and relative mo- 
tion of the two atoms are uncoupled for harmonic traps; 
the weak coupling between them for actual anharmonic 
traps that are used for atom trapping is generally small 
and can be neglected [25^. The potential U bg (R) in Eq. 
0, which describes the relative motion of the atoms, is 
replaced by 



U bg (R) + U trap (R) , 



(15) 



which includes the trapping potential. The continuous 
spectrum of collision energies is replaced by a discrete 
spectrum with trap energies e v and trap cigenfunctions 
4> V (R). Here we take v = 0, 1, ... to represent the eigen- 
values e v > of trap-like states, and v = —1,-2,... to 
represent the bound dimer states with e v < 0. These 
dimer levels of the U bg potential exist even when the trap 
is not present. The continuum function F\ in Eq. is 
replaced by a discrete basis expansion over trap states, 
with time-dependent coefficients since we will be con- 
sidering time-dependent Hamiltonians. The discrete CI 
wavefunction is: 

# n (t) = \l)J2MR)C v (t) + \n)ct> n (R)A n (t). (16) 

v 

Although this equation is written using the " bare" or di- 
abatic basis, it is also possible to diagonalize the Hamil- 
tonian matrix and produce a " dressed" or adiabatic basis 
in which the wave function can be expanded. We will set 
up the coupled equations in the diabatic basis, since that 
is very simple. Below we will show populations calculated 
in both the diabatic and the adiabatic basis sets. 

We assume we start at time to with all the population 
in the ground state of the trap, v = 0, with the magnetic 



4 



field Bo chosen so that e n (Bo) > is above threshold. 
Then the magnetic field is changed linearly in time so 
that the resonance level is ramped below threshold: 



B(t) 



Bq t < t 

if(t-to) t>t Q 



(17) 



The sign of dB/dt depends on the sign of s„ ; clearly 
s n dB/dt < for downward ramps. We only use exam- 
ples in this paper for which s n > and dB/dt < 0, 
and we refer to \dB/dt\ as the speed of the ramp. In- 
serting the CI wavcfunction i|16fl in the time-dependent 
Schrodinger equation gives the coupled set of equations 
that govern the evolution of the CI coefficients: 



ihA r 



ihC v 



[e n (B ) 



8B 

1 at 



t]A n + ^ ' Vn,yC v 



~~\~Vn,vA n . 



(18) 



It is a straightforward procedure to solve Eqs. (|18|) once 
we know the coupling matrix elements 



V n , v = {<Pn\W n ^{R)\^ v ) 



(19) 



These can be evaluated numerically for bound levels be- 
low threshold with v < 0, whereas for trap levels with 
v > they can be found readily from the free-space ma- 
trix elements 



(^ n \W nA (R)\<p v )J— = (<f>„\W n<1 (R)\(t> e ) (20) 
V oe v 

where dv/de v measures the density of trapped states in 
the vicinity of e v e„. Using Eq. i|10[l . we find [l2T| 



r„(e„) de\ 1/2 



2ir dv 



(21) 



There are two cases we will consider. For a har- 
monic potential with frequency Ljf ll and harmonic os- 
cillator length Lh = y^J/MjJh, the density of states is 
a constant, de/dv = 2huif l . If Lh is much larger than 
Ab g , the small energy shifts due to atom interactions j2^| 
can be ignored, and the spherical trap energy levels are 
e v = hujh^ + 2v). The other case is a spherical box of 
length Lb. The box eigenvalues are e v = huJb{v + l) 2 , 
with hwb = (h 2 /2/i)(tt/ Lb) 2 , and the density of states is 
de/dv = 2hL0 b (v + 1). 

Using these properties and Eq. (|13J) . we find for the 
harmonic and box cases: 



huj h ^2V3/7r^/\a h S h \(l + 4w/3) 1/4 , (22) 



V* = hu b V2y/\a b S b \(l + v) , 



(23) 



where a.; = Ab g /Li and <5,; = A n s n /(huji) for i = h,b 
are dimensionless coupling parameters. If the last bound 
state in the potential is not too close to threshold, that 



is, if the scattering length is not too large in magnitude, 
the V n ^\ coupling to the last bound state v = —1 can 
be ignored in Eq. l(T%|) . Thus, a knowledge of the reso- 
nance parameters Ab g , A n and s n is sufficient to calculate 
the time dependent dynamics for a given trap using Eqs. 
CHJ- 

Mies et al. [I3 showed that the populations of the 
ground trap level po and the resonant level p n after a 
linear ramp of resonance energy across threshold are ac- 
curately described by a Landau-Zener (LZ) curve cross- 
ing formula. Assuming an initial population of unity in 
the ground trap level at t = to, the final populations are: 



Po = e 
Pn = 1 



where 



c)B I 
1 dt I 



Li 



Abg A n 



OB 

dt 



(24) 
(25) 



(26) 



and fi = 1 if i = b and fa = V3/tt if i = h. Note 
that the LZ parameter depends on the resonance only 
through the ratio Ab g A n / ^ and does not depend on s n . 
The LZ populations are only asymptotic ones, for long 
time t after the crossing. It is necessary to integrate the 
coupled Eqs. (|18fl in order to obtain the time-dependent 
populations. 

The number of atoms lost from the initially populated 
V = trap ground state level does not depend on the 
direction of the ramp, that is, whether it starts out above 
or below threshold. If the ramp starts with the resonance 
level above threshold and moves it below, then the loss of 
atoms from the v = level represents formation of stable 
molecules in the resonance level. If the ramp starts with 
the resonance level below threshold and moves it above 
into the scattering continuum, then the loss of atoms 
from the v = trap level results in heating, that is, 
population in levels with v > 0. 

As a specific example, Fig. shows the time depen- 
dent populations calculated by integrating the coupled CI 
Eqs. I|18|) for the model shown in Fig.^ As the resonance 
nears the crossing at time t = 0, the populations undergo 
transient oscillations and settle down at late time at the 
values given by the Landau-Zener model. If we apply the 
Demkov-Osherov |26l [23] linear multiple curve crossing 
model to this case, for a ramp starting at t = —00, there 
is exactly zero final population transfered to excited trap 
levels, i.e, those with v > 1; such populations can only be 
excited transiently around the time of the crossing. All 
final population either remains in v = or is transfered 
to the resonance state n. Ref. ^3] showed that for either 
downward or upward ramps, the transition probablilites 
could be calculated accurately for one or for a sequence 
of curve crossings using the Landau-Zener expressions. 
This is because for a linear set of curve crossings, the 
Demkov-Osherov model is an exact representation of the 
quantum dynamics. 
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FIG. 2: Time dependent populations projected in the dia- 
batic and adiabatic basis sets for the model problem shown 
in Figure for a ramp speed of 0.5 mT/ms. The ramp 
starts at to = — 25oj^~ 1 . The diabatic levels shown are the 
"bare" ground state v — (2 atoms) and the resonance level 
n ("bare" molecule). The adiabatic populations shown are for 
the lowest two adiabatic states in Fig. Initially, the low- 
est adiabatic level is almost indistinguishable from the "bare" 
ground trap level. This level adiabatically evolves to become 
nearly indistinguishable from the "bare" resonance level for 
times long after the threshold crossing at t = 0. The initially 
unpopulated v = 1 trap level adiabatically evolves to become 
the v — ground trap level long after the crossing. Although 
the final state populations evolve to the Landau Zener limit, 
the transient populations near t = may in the general case 
depend on the size of the CI basis set, especially in cases where 
the background scattering length is large in magnitude. The 
populations strongly depend on the basis in which they are 
projected. 



populations at a frequency determined by the binding 
energy of the dressed molecular resonance state. Ref- 
erence [3(j applied a model of two atoms interacting in 
a trap to the condensate experiments of Ref. They 
adjusted their model parameters to account for the qual- 
itative features of the experiment. The next section de- 
scribes the alternative approach introduced by Ref. [IJ 
to incorporate the many-body physics of a condensate 
into a model of trapped atom pairs. This model has no 
adjustable parameters. 

87 Rb 100 kHz trap, 6-level model 



,, , f\ c m i -- Diabatic basis 

Ramp speed = 0.5 T/s _ Adiabatic basis 




FIG. 3: Energy levels for the same 6-level model described in 
Fig. Q except with a backwards ramp of the resonance level 
above threshold at the same ramp speed after a hold time 
St at a constant magnetic field, where the diabatic resonance 
energy is held at an energy of — lOhuih- If the ground trap 
state v = is initially populated, then final populations will 
be distributed among the trap final states with v > 0. 



A variety of different ramp types can be applied. Fig- 
ure |3 gives an example of a double ramp, which first 
ramps the resonance level across threshold, then after 
a hold time St, ramps it back in the opposite direction 
to be far above threshold at long time. Figure 0| shows 
how the final adiabatic populations at late time oscillate 
as St is changed. The populations in the ground v = 
and first excited v = 1 trap levels oscillate with a period 
t = 2-k/lo given by the difference hcu between the energies 
of the (dressed) adiabatic trap ground level and molecu- 
lar resonance level during the hold interval. For example, 
measuring the oscillation frequency for atom pairs in op- 
tical lattice cells for different hold magnetic fields could 
be used to measure s n = de n /dB. 

The oscillations in Fig. 0] for the case of two atoms in 
a lattice cell bear a similarity to observations on a 85 Rb 
BEC using pulsed magnetic field ramps with a different 
time dependent configuration Q, Il4j . In these experi- 
ments the adiabatic bound state approached threshold 
from below and remained below threshold at all times. 
However, the timc-dcpcndcnt ramps induced nonadia- 
batic transitions that resulted in the formation of dressed 
molecules and oscillations in the final atomic condensate 
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FIG. 4: Final state populations at long time t for the 

ramp sequence shown in Figure^ Populations are shown for 
the adiabatic levels with the three lowest energies in Figure 01 
The oscillation frequency is determined by the energy differ- 
ence between the two lowest adiabatic levels (the approximate 
molecular resonance state and the trap ground state) during 
the hold time St. 
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IV. FESHBACH RAMPS IN A BEC 

Mies et al. 01 devised an approximate procedure to 
take into account ramping a time-dependent resonance 
level across threshold in a BEC. They used the theory 
of atom pairs in the trap developed in Section IIIII and 
showed what changes must be made to account for the 
many-body effects in the condensate. There are two es- 
sential effects that must be treated. First, assuming that 
the background scattering length is not close to zero, the 
many-body mean field interactions cause the condensate 
size to expand to occupy a much larger volume than a sin- 
gle atom in the harmonic trapping potential, which would 
be confined to a volume of radius equal to V3 times the 
harmonic length for a single atom, L a = L^/^/2. In such 
a case, the condensate is characterized by the Thomas- 
Fermi radius 



where 



Ltf = L a ?ft , 



(27) 



(28) 



Since the decay width in Eq. I|10|) . which is needed in 
Eq. I|21[) . is proportional to e 1 / 2 , we need a criterion for 
picking a value of collision energy e app ropriate for the 
ground state of a condensate. Ref.[l2| used the result 
that the mean kinetic energy per atom in the condensate 

is HEl 



7? 



huj h 



(29) 



where uJh is the mean harmonic frequency of the confining 
potential. Since the mean field potential exactly cancels 
the harmonic potential in the Thomas-Fermi approxima- 
tion to make a flat effective potential for distances less 
than Ltf, we use for the ground state of the relative mo- 
tion of a pair of condensate atoms the ground eigenstate 
of a box that has a zero point energy huj^ equal to (T a ) 
in Eq. (|29|1 (the kinetic energy of relative motion is equal 
to that of a single atom of the pair) . This gives 



L h = L 



TF 



ln(- 



(30) 



2683 1 



Therefore, the harmonic trapping potential for the atom 
pair is replaced by a box potential that has the size Lb. 

The second effect of having a condensate is that the 
ground state is macroscopically occupied by N atoms. 
Each atom has N — 1 atoms with which it can interact. 
Since N 1, Ref. 0] introduced the factor y/N into the 
V n fi matrix element in Eq. (|23|l . The resulting matrix 
element is proportional to \J N /L\. Consequently, the 
Landau-Zener parameter in Eq. I|26|) is proportional to 
an atomic density. 



If we assume a T = condensate, there is no initial 
population in the excited trap levels with v > 0. Con- 
sequently, the dynamics is dominated by coupling of the 
resonance to the ground v = trap level, which is en- 
hanced by the ^/N factor, and we can ignore the excited 
levels, unless they become macroscopically populated in 
the course of the dynamics. As discussed in Section llVl 
there are no transitions from v = to v > levels by 
a downward ramp, since these transitions are in a coun- 
terintuitive direction and are rigorously forbidden in a 
Landau-Zener model. Thus, we might expect that the 
excited levels can be neglected for a downward ramp. 
However, they can be populated when there is an up- 
ward ramp that moves population from the ground state 
to higher energv. This is what ha ppe ns in the MIT ex- 
periments 0, modeled in Ref. |l2l | . Even in this case 
the removal rate of atoms from the ground state depends 
only on V n fi and does not depend on the excited states. 
We will show in the next Section that the predictions 
of this simple picture compare well with the results of 
proper many-body calculations of the system. 

To calculate the time-dependent loss of condensate 
atoms when the Feshbach resonance is ramped across 
threshold from above or below, given that trap excita- 
tion can be neglected, all that is needed is to use the 
two coupled Eqs. (|18|l for the resonance state and the 
V = trap state, with the modified V n $ matrix ele- 
ment of Eq. (|23[) multiplied by \f~N and using the box 
length from Eq. I|3U|) . These coupled equations describe 
an atom pair in the condensate. If N is the total num- 
ber of atoms in the system, the number of atom pairs 
is N/2. Of these atom pairs, after a downward ramp 
the fraction po represents BEC atom pairs and the frac- 
tion p n = 1 — po represents molecules. Consequently, the 
number of remaining condensate atoms is Npo and the 
number of molecules is 7V p n / 2, with Np n atoms com- 
bined into molecules. Ref. [Tj] showed that the Landau- 
Zener formula gives an excellent account of the transition 
probabilities. Let us call this method the LZ method. 

It is straightforward to make a relatively minor cor- 
rection to the theory to account for the nonlinear na- 
ture of the dynamics. In the case of a BEC, we found 
that it is necessary to multiply the matrix clement V£ 
in Eq. I|23f) by y/~N. However, N decreases dynamically 
during the course of the ramp as atoms are transferred 
out of the condensate. This is readily treated by intro- 
ducing the time dependent number of condensate atoms, 
N\C (t)\ 2 m Eqs. (fT5fl for the condensate case, including 
only the v = ground state of the trap: 



ifiA n 
ihCo 



e C + V n , ^/N\CtfA n . (31) 



This makes the dynamics nonlinear, as in the Gross- 
Pitaevsky equation, and corrects the LZ model in the 
right direction. Let us call this the CI-BEC method. 
Given the prescription in Eq. (|30|) , based on a physically 
motivated choice of collision energy, these equations have 
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no adjustable parameters. 

If the ramp is upward, the resonance energy ramps 
from below threshold to higher energies. Even if the cou- 
pling is very strong between the ground trap (i. e., con- 
densate) level and the resonance level, the coupling will 
be much weaker to the individual v > 1 levels because 
they are initially unoccupied. So the crossings with the 
higher trap levels will have very small coupling and tend 
to be diabatic. Thus there is little loss to the first few 
excited trap levels, and the population is carried to very 
high levels before it finally decays into the quasicontin- 
uum of trap levels to produce hot atoms. For the special 
case when only molecules are initially present without 
condensate atoms, a very simple expression 0, Q f° r 
the energy distribution of hot atoms can be worked out 
by simply assuming Fermi-golden rule decay of the reso- 
nance level n at a time-dependent rate T n (e n (t))/h given 
by Eq. HSJ. 



V. SPECIFIC BEC CALCULATIONS 

This Section will summarize the results of some specific 
calculations involving Fcshbach ramps for 23 Na, 87 Rb, 
and 133 Cs condensates. 



A. 23 Na F = 1, M = +1 

Figure El shows LZ and CI-BEC calculations for the 
condensate atom loss 1 — po, using the 85.2 mT and 90.7 
mT resonances for collisions of 23 Na F = 1,M = +1 
atoms. Mies et al. H2 had previously applied the LZ 
method to these resonances to explain qualitatively the 
upward ramp experiments by the MIT group, which had 
measured the loss of condensate atoms associated with 
different ramp speeds dB/dt (13 . The three model 
parameters arc A^ g = 3.32 nm, A„ = 0.1 mT, and 
s n /h = 52.9 MHz/mT for the broad 90.7 mT resonance 
and A{,g = 3.38 nm, A„ = 1 /iT or 0.25 /xT, with the 
same s n for the narrow 85.3 mT resonance. The cal- 
culated width and estimated experimental width agree 
for the former, but differ by a factor of 4 for the lat- 
ter 1^3. The calculations assumed a trap with the ex- 
perimental mean trap frequency of (Jh = 27r(700 Hz) and 
TV = 900000 atoms. The actual experimental optical trap 
was not spherical, but cigar-shaped. 

Figure [5] shows that the model gives the right order 
of magnitude of the atom losses. If the renormalization 
of the matrix clement for many-body effects described 
in Section llVl were not taken into account, and a "bare" 
harmonic trap were used, then the ramp speeds would 
differ by about three orders of magnitude to give the same 
loss fractions. The CI-BEC calculations agree with the 
LZ calculations for fast ramp speeds (small inverse ramp 
speeds) where the atom losses are small, but give lower 
atom losses than the LZ calculations at slow ramp speeds, 
where the atom losses are large. One would expect the 



LZ model to be better for fast ramps with small atom 
losses, since the initial condensate is less disturbed by the 
process. The CI-BEC calculations for the fastest ramp 
speeds are closer to the experimental losses for the narrow 
85.3 mT resonance when the theoretical width 1 fiT is 
used instead of the estimated experimental width 0.25 
of Ref. |3. 

Thus, it seems that the LZ formula in Eq. [25] is useful 
for making order-of-magnitude estimates of atom losses 
using only three experimentally accessible parameters, 
Ab g , A n and dB/dt, and two trap parameters, <Jh and 
N, using the prescription in Eq. (|30|l for the effective box 
size. It is clear that proper many-body calculations are 
needed for quantitative studies. 
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FIG. 5: Atom loss versus inverse ramp speed \dB/dt\~ l for 
upward ramps of the 85.2 mT and 90.7 mT resonances in 
23 Na F = 1,M = +1 collisions. The figure compares the 
calculated LZ (dashed line) and CI-BEC (solid line) results 
with the experimental data points I32I . In the case of the 
narrower 85.2 mT resonance, the CI-BEC results are shown 
for two different values of the resonance width, A„ = 1 fiT 
and 0.25 f/I. 



B. 



^ Rb F = 1, M = +1 



Goral et al. ^3 carry out detailed calculations for the 
100.76 mT resonance in 87 Rb F = 1, Al = +1 collisions, 
comparing many-body calculations based on the method 
of Kohler et al. [3 El with the LZ and CI-BEC models 
described here. This is the same resonance discussed in 
Figs. 1 1141 for two atoms in an optical lattice cell. Our 
many-body methods permit a clear way of counting the 
populations of the condensed and non-condensed atoms 
and the atoms combined in a molecular bound state. The 
latter refers to the adiabatic, or "dressed" , bound state of 
the two-body Hamiltonian, not the diabatic, or "bare", 
bound state that is only an intermediate theoretical con- 
struct and not an exact two-body eigenstate. Thus, the 
method accurately counts the number of molecules that 
could be measured in an experiment in which they are 



captured. We will summarize some of the results of 
Ref. here. The key finding is that the CI-BEC model 
gives remarkably good results when compared to many- 
body calculations. 

Fig[f)lgives an example of calculations for a uJh = 27r(10 
Hz) trap with 50000 atoms. The ramp starts with the res- 
onance detuned far above threshold with all population in 
the condensate ground state and ends with the resonance 
tuned far below threshold. The LZ and CI-BEC methods 
predict the number of remaining condensate atoms (Npo) 
and number of atoms combined in molecules (Np n ) in 
good agreement with the many-body calculations |34j . 
In this case, the fractional depletion of atoms is modest, 
and the LZ and CI-BEC methods only differ slightly in 
their predictions. Figure \7\ gives an example of a ramp 
of the same resonance in a uJh = 27r(100 Hz) trap. There 
is much more atom depletion and molecule formation in 
the 100 Hz trap than the 10 Hz trap for the same ramp 
speed, because the tighter trap gives larger matrix ele- 
ments between the resonance state and the trap ground 
state. The CI-BEC method agrees well with the many- 
body calculations, whereas the LZ method predicts too 
much atom loss from the condensate for the slower ramp 
speeds. Thus, the very simple LZ model is good for es- 
timating the magnitude of the ramp speed needed for 
conversion of atoms to molecules, although the actual 
efficiency of molecule production saturates less rapidly 
tha 




Inverse ramp speed (s/T) 



FIG. 6: Final populations versus inverse ramp speed 
\dBfdt\~ 1 for a downward ramp of the 100.76 mT resonances 
in 87 Rb F — 1, M = +1 collisions in a 10 Hz trap with 
N = 50000 atoms. The figure compares the calculated LZ 
(dashed line) and CI-BEC (solid line) results with the many- 
body calculation for the remaining condensate atoms (solid 
circles) and atoms combined in the molecular resonance state 
(solid squares). 



C. Cs F = 3, M = +3 

Herbig et al. @ have succeeded in using a Feshbach 
ramp to produce Cs2 dimer molecules from a 133 Cs con- 
densate of F = 3, M = +3 atoms trapped in an optical 




Inverse ramp speed (s/T) 



FIG. 7: Final populations versus inverse ramp speed 
\dBfdt\~ 1 for a downward ramp of the 100.76 mT resonances 
in S7 Rb F = 1,M = +1 collisions in a 100 Hz trap with 
N = 50000 atoms. The figure compares the calculated LZ 
(dashed line) and CI-BEC (solid line) results with the many- 
body calculation for the remaining condensate atoms (solid 
circles) and atoms combined in the molecular resonance state 
(solid squares). 

trap [33. We have calculated the properties of the 2 
mT resonance in these experiments using the theoreti- 
cal model of threshold scattering developed by Leo et 
al. 0| calibrated from the Feshbach spectroscopy exper- 
iments of Vuletic et al. . We have also calculated the 
molecule production rates using the LZ, CI-BEC, and 
many body methods, using our calculated resonance pa- 
rameters, Ab g = 8.63 nm, A„ = 0.5 [iT, and s n /h = 7.98 
MHz/mT. The s n value is in excellent agreement with 
the measured magnetic moment of the molecules Q . We 
will explain our calculations in full details elsewhere. 

The diabatic resonance level \n) has (/to/ ,lmi) = 
(44, 42) molecular symmetry, where / and I respectively 
represent spin and relative angular momentum with pro- 
jections mf and to/. This g-wavc symmetry resonance 
state \n) is only very weakly coupled to the entrance 
channel s-wave collision. Our calculated resonance width 
is very sensitive to the short-range second-order spin 
orbit corrections to the effective spin-spin interaction, 
and is subject to error, possibly as much as a factor of 
two. We find fairly good agreement in the limit of fast 
ramp speed with preliminary data from the Innsbruck 
group [36j on the number of molecules produced versus 
inverse ramp speed for a cJ/, = 27r(ll Hz) trap with 50000 
atoms. 

There is also a 133 Cs resonance for F = 3, M = +3 
collisions near 4.8 mT with (frrif,lmi) = (44,22) molec- 
ular symmetry. Figure [5] shows our calculated scattering 
length between and 10 mT for collisions of atoms in this 
state, using the model of Leo et al. |l7j including s- and 
d-basis functions in the wavefunction expansion. Atom 
losses due to ramping across this resonance were briefly 
described in the original condensate experiment |35j , but 
with insufficient information for us to model. It should 
also be possible to make molecules with this resonance. 
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The resonance parameters are A\, g = 47.9 nm, A n = 15 
fiT, s n /h = 20.9 MHz/mT. 



133 Cs Scattering length, F=3, M=+3 
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FIG. 8: Scattering length versus B, calculated using i = and 
2 basis functions, with the inset showing the (44,22) resonance 
at 4.8 mT. The (44,42) resonance at 2 mT studied by Herbig et 
al. is missing, since i = 4 basis functions were not included 
in this calculation. 
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FIG. 9: Diabatic and adiabatic energy level versus B — B„ 
for a CI model of a pair of Cs F = 3, M = +3 atoms in 
a 529 nm box. The points show the exact numerical eigen- 
values using the full spin Hamiltonian in the model of Leo et 
al. [TtI . The arrow shows the predicted location of the free- 
space threshold crossing, using the shift formula Eq. 11411 . 

The threshold crossing of the 4.8 mT resonance is sig- 
nificantly influenced by coupling to the last bound state 
of the background potential, which has a energy of only 
t-i/h — —45 kHz. Figure [5] shows the bound states of 
relative motion of two F = 3, M = +3 133 Cs atoms in a 
529 nm box (529 nm = 10 4 atomic units). The lines show 
the diabatic and adiabatic energy levels calculated using 
the CI model described in Section ITTll The figure shows 
these are in good agreement with the numerical eigen- 
vaulcs calculated from the full spin-dependent Hamil- 
tonian and molecular potentials of the Leo et al. ^3] 



model. The coupling of the resonance level |n) to the 
last bound state \v = —1) of the backgound potential 
is also included, using the calculated matrix element 
V n ,-\/h = 183 kHz. The strong avoided crossing of the 
|n) and \v = —1) levels is evident in the figure. This 
avoided crossing is responsible for the large shift in the 
actual threshold crossing position. The shift for this trap 
case is close to the free-space shift of — 1.10A„ predicted 
by Eq. i|14|) for our model parameters, Af, g = 47.9 nm 
and A = 5.11 nm. 

Figure shows the adiabatic energy levels for the LZ 
Hamiltonian matrix using a two-level model with |n), 
|u = 0) in the basis and a three-level model with \v = —1) 
also in the basis. In the two-level model there is an 
avoided crossing centered on the time that the resonance 
level \n) crosses threshold, and adiabatic evolution takes 
the initially populated trap ground state atoms in level 
\v = 0} to molecular resonance level |n) population for 
times long after the crossing. By contrast, the strong 
interaction with the last bound state in the three-level 
model shifts the threshold crossing to later time, and 
adiabatic evolution would lead to producing molecules in 
the \v = —1) molecular state instead of in the diabatic 
resonance state \n). 

Figure 1111 shows the time-dependent populations cal- 
culated using the CI-BEC model for the case shown in 
Fig. Since the LZ parameter describing the loss only 
depends on the V u q matrix element, the LZ model pre- 
dicts the same final loss of initial atomic population for 
both the two- and three-level models. However, the de- 
tailed population dynamics, and the nature of the fi- 
nal molecules produced, are profoundly different between 
the two models. When we consider the effect of the 
last bound state, the actual transfer of population from 
atoms to molecules occurs near the time of actual cross- 
ing, shifted by the strong interaction with the last bound 
state, and results in the formation of molecules in the 
last bound state \v = —1) of the background potential, 
not formation of resonance state molecules It is in- 
teresting to note that in the three-level model, even at 
intermediate times near the resonance crossing, there is 
never significant population in the resonance state \n). 

In closing, we note that experiments could also be done 
for 133 Cs F = 3, M = +3 atom pairs trapped in op- 
tical lattice cells, as described in Section IIIII Similar 
qualitative behavior to that described in Figs. EH an d 
1111 should occur. The dynamics would contrast sharply 
from that described in relation to Figs. ^ and |21 for 87 Rb 
F = 1,M = +1, for which At, g is much closer to A, 
and the last bound state is much more deeply bound by 
e-i/h = -23 MHz [l6j. Interaction of |n) with \v = -1) 
results in only a slight modification of the threshold cross- 
ing, and would not significantly change the results in 
Fig. |21 for ramps that stop before allowing e n to cross 
e_i. Of course, if the ramp should carry e n to be close 
to or to cross e_i, the mixing of these levels leads to 
the very interesting physics that has been exploited by 
Diirr et al. £| in making molecules from 87 Rb condensate 
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atoms. 



VI. SUMMARY 
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FIG. 10: Energy levels of the LZ Hamiltonian for a XS6 Q,s, 
F = 3, M = +3 condensate in a uJh = 27r(20 Hz) trap with 
12500 atoms, for the cases of two-level and three-level models 
for a downwards ramp near 4.8 mT with speed 0.4 T/s. The 
dashed lines show the diabatic levels. The numbers label the 
adiabatic levels in order of increasing energy. The three-level 
model also includes the last bound level \v = — 1) of the open 
channel potential at -45 kHz. The strong interaction between 
the resonance level \n) and \v = —1) leads to a strong avoided 
crossing that changes the nature of the threshold crossing, 
shifting it to 88 /is later in time for the ramp speed used here. 
In the three-level model the adiabatic level 2 evolves from its 
initial character as the trap ground state \v = 0) to its final 
character as the \v = — 1) molecular state. 
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FIG. 11: Calculated diabatic populations using the CI-BEC 
equations for the two- and three- level models shown in Fig. 1101 
The inclusion of the last bound state profoundly changes the 
detailed dynamical evolution, although the final state popula- 
tions of condensate atoms and molecules is the same in both 
cases and is close to the value given by the LZ formula. In 
the three-level case, the population of the resonance state |n) 
always remains small, and is barely visible in the figure. 



We have presented a review and two extensions of the 
Feshbach ramp model of Mies et al. 0| for describing 
the conversion of pairs to trapped atoms to molecules by 
ramping the energy of a Feshbach resonance state from 
above to below a collision threshold. The review out- 
lines the key ideas behind the method. The extensions 
arc given by in Eqs. ij^ - ij^ and lj!Tj"|) . The former 
equations express the matrix elements using the Wigner- 
law form of the coupling, and the latter Eq. Q31[) intro- 
duces nonlinear dynamics for the case of condensate evo- 
lution. Given the number of atoms in the trap and the 
mean trap frequency, the model is fully characterized in 
terms of only three nonadjustable resonance scattering 
parameters which can be measured or calculated: the 
background scattering length, the width and the mag- 
netic moment of the resonance state. The model can be 
supplemented by adding the coupling matrix element to 
the last bound state in the background potential if nec- 
essary. Given the ramp speed, the model predicts the 
loss of atoms from the initially populated ground state 
of the trap and the conversion efficiency to molecules. 
Simple estimates based on the Landau-Zener curve cross- 
ing formula give the correct order of magnitude of atom 
loss when compared to experimental cases. The CI-BEC 
equations compare quite favorably with the predications 
of many-body theory for condensates. 

The model can be extended to other cases. It could be 
readily adapted to one- or two-color photoassociation. 
This can be treated by a Feshbach resonance formalism 
that includes the decay of the excited state by spon- 
taneous emission. In the two-color case, it is possible 
to su ppr ess this decay by using an interference mecha- 
nism [3jfl , The model also applies to two Fermions with 
different spins, and could be applied without modifica- 
tion to pairs of fermionic atoms in single trapping cells. 
It is possible that the model could be adapted to quan- 
tum degenerate mixed fermionic gases. One important 
physical effect that the model leaves out is the role of 
other collisions during the ramp process. For example, 
if a third atom approached the colliding pair of atoms, 
a resonance-modified three-body process occurs. There 
clearly is a rich variety of phenomena that can be stud- 
ied based on manipulating resonant scattering processes 
in optical lattices or quantum degenerate gases. It is 
hoped that models like the one presented here will help 
stimulate thinking about these phenomena. 
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